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Abstract 

The rovibrational kinetic energy for an arbitrary number of rigid molecules is com- 
puted. The result has the same general form as the kinetic energy in the molecular 
rovibrational Hamiltonian, although certain quantities are augmented to account for 
the rotational energy of the monomers. No specific choices of internal coordinates or 
body frame are made in order to accommodate the large variety of such conventions. 
However, special attention is paid to how key quantities transform when these conven- 
tions are changed. An example system is explicitly analysed as an illustration of the 
formalism. 

1 Introduction 

The rovibrational kinetic energy of a molecular (rigid body) complex has previously been 
computed explicitly for a number of specific cases. Examples include calculations for the 
atom-monomer system, consisting of a point particle and a rigid body, by Brocks and van 
KoevenJlJ, van der Avoird ||, and Makarewicz and Bauder 0; the molecular dimer, con- 
sisting of two rigid bodies, by Brocks et al. [[| and van der Avoird [[|; the molecular trimer 
by Xantheas and Sutcliffe|| and van der Avoird, Olthof, and Wormer [7|; as well as closely 
related systems of single molecules with internal rotation M, R|. In the present article we 
derive the rovibrational kinetic energy of an arbitrary molecular complex containing an ar- 
bitrary number of rigid bodies. That is, we express the kinetic energy in terms of the total 
angular momentum of the complex and the momenta associated with the internal degrees 
of freedom. The term "rovibrational" is perhaps misleading since it implies small amplitude 
vibrations of the complex, an assumption we do not make. A rovibrational kinetic energy 
for a general molecular complex has been proposed earlier by Makarewicz and Bauder [[|; 
their analysis differs markedly from ours in that they do not impose rigidity conditions on 
the monomers nor do they distinguish the relative rotational motion of the monomers from 
the internal vibrations of the monomers. 

One method for deriving the rovibrational kinetic energy of a system of rigid bodies 
would be to begin with the rovibrational kinetic energy of a system of point particles and 
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then impose rigidity constraints within certain subsets of these particles. A general analysis 
of internal constraints on n body systems has been given by Menou and Chapuisat [TIJ and 



by Gatti et al. ||11|| . These authors observe that such a formalism may be applied to find 
the rovibrational kinetic energy of a rigid body complex, though they do not derive such a 
kinetic energy. The constrained systems approach has the advantage of naturally allowing 
the relaxation of the rigidity constraints to include small internal vibrations of the monomers. 
However, within a strictly rigid body, the positions, masses, and velocities of the constituent 
particles are irrelevant. Rather, only the overall orientation, moment of inertia, and angular 
velocity are of interest. Therefore, in our derivation of the rovibrational kinetic energy, 
we assume from the outset that the bodies are rigid and use a kinetic energy consisting of 
only the translational and rotational energies of the rigid bodies. It should be mentioned 
that whichever approach is followed, there exists an ambiguous extrapotential term in the 
quantum kinetic energy. The origin of this ambiguity rests in the lack of knowledge about the 
potential which confines the system to the manifold of rigid shapes. See for example Kaplan, 
Maitra, and Heller [ I2"| . For simplicity, we neglect any extrapotential terms arising from the 
constraint process and adopt the standard form (translational plus rotational energies) for 
the quantum kinetic energy of a rigid body. 

Our approach is modeled on the derivation of Ref. valid for clusters of point particles, 
but is augmented to include the rotational kinetic energy of the monomers; the current paper 
therefore generalises many results valid for point particles to systems of extended rigid bodies. 
Also, in deriving the rovibrational kinetic energy, we make no specific choice of internal 
coordinates or body frame. This allows our results to be applicable to the wide range of 
coordinate and frame conventions suitable to different molecular complexes. We discuss how 
the various quantities appearing in the kinetic energy transform under changes of internal 
coordinates and body frame. We also present two distinct decompositions of the kinetic 
energy into rotational and vibrational (internal) contributions. The first decomposition 
is independent of the conventions for internal coordinates and body frame. The second 
decomposition has a form more common to rovibrational Hamiltonians in the literature, 
though it is not independent of coordinate and frame conventions. 

The current paper has been influenced by recent work exhibiting the importance of 



gauge theory and geometric phase in the study of rovibrational coupling O, [TJ, ITR [16 



17| , [18], [1{J. In addition to rovibrational coupling, the geometric, or Berry's, phase |£0| has 
important applications in Born-Oppenheimer theory [Ej]], opticsf^, |2"3|] , and guiding centre 
motion p4j, to name but a few examples. The influence of gauge theory on our analysis 



here is most readily evident in the form of the kinetic energy Eq. (0), which owing to the 
appearance of a gauge potential is reminiscent of the kinetic energy of a charged particle in a 
magnetic field. Notions of gauge invariance and covariance have also influenced our analysis 
of transformation properties presented in Sect. [2.3[ Although our development is influenced 
by the techniques and concepts of gauge theory and geometric phase, this paper requires no 
specific background in either field. 

The article is structured as follows. Section |2| is devoted to the classical kinetic energy. 
The principal computation of the article, culminating in the rovibrational kinetic energy 



Eq. (|27|), is contained in Sect. 2.1. In Sect. 2.2 we show how the rovibrational kinetic energy 



of Eq. ( P7| ) can be placed in a different form, that of Eq. fl3"3"D, which is more common in 
the rovibrational literature. Sect. 0|is a systematic discussion of how important quantities 



2 



transform under changes in the internal coordinates and body frame. Sect. |3| is an aside in 
which we discuss the form of quantum kinetic energy operators in general. In particular, 
Sect. |3J] discusses an alternative to the usual Podolsky form of the quantum kinetic energy, 
and Sect. |3.2| illustrates how scaling the wave function introduces an extrapotential term into 
the kinetic energy. In Sect. |] we apply the results of Sect. |3| to compute the quantum kinetic 
energy of a molecular complex. We work through a concrete example in Sect. || applying our 
formalism to the system of a single monomer and an atom. Sect. ^ contains conclusions. 



2 The classical kinetic energy of a molecular complex 
2.1 The principal derivation 

In this paper, a molecular complex is modeled by a collection of n rigid bodies or monomers. 
No constraints are placed on the positions or orientations of the monomers or on the sym- 
metry of the moment of inertia tensors. In particular, we allow the monomers to be point 
particles (atoms), collinear bodies, or noncollinear bodies. However, for mathematical sim- 
plicity, we initially assume that each monomer is noncollinear. Then, after deriving the 
rovibrational kinetic energy, we comment on the straightforward generalisation of allowing 
complexes containing collinear monomers and point particles. 

We begin by introducing three classes of frames which are important in our derivation. 
The space frame (SF) is the inertial, laboratory-fixed frame. There are n individual body 
frames (IBF), one for each body; IBF a, a — 1, . . . , n, is fixed to rigid body a and rotates 
with the body. Finally, there is a single collective body frame (CBF), which is fixed to the 
complex as a whole. The CBF differs from the other frames in that it must be specified for 
each shape or internal configuration of the complex. In general, this specification produces 
singularities in the CBF, a fact which has been studied in three- and four-atom systems 
p5| , p6| . There is no canonical way of choosing the CBF, though there are several methods 
commonly used in such Hamiltonians, such as fixing the CBF to the principal axes of the 
complex or to one of the IBFs. In this paper we will make no specific choice of CBF. 

We employ the following notation to denote the frame to which the components of a vector 
are referred. For an arbitrary vector v, an s superscript, that is v s , indicates components in 
the SF; an ia superscript indicates components in the IBF of body a; and a c superscript 
indicates components in the CBF. (For notational simplicity, at a certain point we will drop 
the c superscript, leaving it understood thereafter that a vector without a superscript is 
implicitly in the CBF.) The components of v in the various frames are related by proper 
orthogonal 3x3 matrices, which we define by 



v 

V s ' 

v c 



Rv c , 

- J a v 



(1) 

(2) 
(3) 



The matrix determines the orientation of body a, that is, the IBF of body a, in the SF; 
the matrix R determines the orientation of the entire complex in the SF; and the matrix S^, 
determines the orientation of body a in the CBF. 
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The configuration of the complex of n rigid bodies is fully specified by the centre of mass 
position of each body and the orientation of each body. To eliminate the overall translational 
degrees of freedom, we fix the centre of mass of the entire complex at the origin. The centre 
of mass positions of the n bodies are then determined by n — 1 Jacobi vectors E7I, IT3 



r*, a = 1, ... ,n — 1. The orientations of the rigid bodies are specified by the n matrices 
G £0(3), a = 1, ... ,n. Taken together, r* , a = 1, . . . , n — 1 and S* , a = 1, . . . , n specify 
a lab description of the configuration. To shift to an internal-external, or shape-orientation, 
description of the configuration, we introduce 6n — 6 internal, or shape, coordinates g M , 
fi = 1, . . . , 6n — 6. The g M may be separated into 3n — 6 coordinates parametrising the 
distances between the rigid bodies and 3n coordinates (Euler angles) parametrising the 
orientations of the bodies in the CBF. However, here we allow the g M to be completely 
arbitrary, so long as they are invariant under rotations of the complex. The Jacobi vectors 
r£,, referred to the collective body frame, and the matrices S° are both functions of the 
internal coordinates g M . In practice, one may define the CBF by specifying the functions 
v c a (q) and S c a (q), where q without a superscript refers to the collection of all coordinates 
The orientation R G 5*0(3) (defined in Eq. ([!])) of the complex together with the shape q^ 
fully determine the configuration in the SF, as may be seen by the following equations 



< = Rr£(g), (4) 

S* = RS°(?). (5) 

Equation (f|) is an application of Eq. ([]]), and Eq. (|) is implied by Eqs. dH)-©. The above 
equations relate the internal-external description of a configuration (in terms of R and q) to 
the original description (in terms of r* and S*). 

We assume that the masses of the n bodies have been absorbed into our definition of the 
Jacobi vectors. Thus, the kinetic energy of the complex, without the overall translational 
contribution, is 

n— 1 , n 

T =2El^| 2 + 2E<- M >- ( 6 ) 

a=l a=l 

where the dot is used for time derivatives, is the moment of inertia of body a in the SF, 
and uj s a is the angular velocity of body a. 

In general, an angular velocity is a vector which measures the rotation rate of one frame 
with respect to another frame. The components of this vector may be referred to either of 
these two frames (or to an arbitrary third frame for that matter). For example, the angular 
velocity u;* measures the rotation rate of the IBF of body a with respect to the SF; its 
components are referred to the SF. 

We introduce the notation vx for the antisymmetric matrix which maps a vector u into 
the vector v x u. Then, 



w!?x 



rs cs T 

S sT S s 

-'a J a) 



(7) 
(8) 
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where the T superscript denotes the matrix transpose. Note that these two formulas are 
consistent with the change of basis relation Eq. (g), that is uj s a = S*u>^*, as may be seen 
from the general relation 

Q(vx)Q T = (Qv)x, (9) 

where the vector v is arbitrary and Q 6 50(3). 

We proceed by expressing the kinetic energy Eq. (^) in terms of the internal velocities q^ 
and the total angular velocity u> c . The total angular velocity measures the rotation rate of 
the CBF with respect to the SF. It has components in both the SF and the CBF which are 



u s x = RR T , (10) 
o; c x = R T R. (11) 

These equations are analogous to Eqs. (J?]) and (§). Equation (pTT|) permits the time derivatives 
of Eqs. (|J) and (||) to be expressed as 



= Rr^ + Rr^ = R(o; c x^ + r^), (12) 
S* = RS^ + RS^ = R[Kx)S^ + S^], (13) 

where the ", /z" subscript denotes the derivative with respect to the coordinate q^. In the 
above equations, we have used the convention, which we adopt for the remainder of the 
paper, that the Greek indices /i, z/, . . . are implicitly summed from 1 to 6n — 6 when repeated. 
However, the Greek indices a, f3, . . . which label either the Jacobi vectors or monomers are 
summed explicitly. Combining Eqs. (0) and (p~3f ) and using Eq. (|), we find 



<x = R[Kx)+S^V]R T . (14) 
Since S^S^ 7 " is antisymmetric, we define a vector T afl such that 



T a^ x — S^S^ T , (15) 
r%x = Sf S^. (16) 

By comparison with Eqs. (0) and (||), we see that q^T c a ^ is the angular velocity of the IBF 
a with respect to the CBF. By inserting Eq. ([TJ|) into Eq. flT4] ) and using Eq. ([5]), we find 

< = R(^+^;)^ s + ?' 1 t;. (17) 

The above equation expresses the angular velocity of the IBF of body a with respect to the 
SF as the sum of the angular velocity of the CBF with respect to the SF plus the angular 
velocity of the IBF of body a with respect to the CBF. However, this decomposition has no 
inherent physical meaning since it depends on the convention used to define the CBF. By 
appropriately changing this convention, either of these two terms could be made to vanish. 
Using Eqs. ( |I2"D and ( |T7| ) in Eq. (|5]), we write the kinetic energy as 
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T = ~u> • Mo; + a^ • w(f + hi^q\ (18) 



where 



n— 1 n n— 1 



a=l a=l a=l a=l 

n— 1 n n— 1 n 

J2 TaX r ^ + MqT ^ = r « x r ^ + SqM « ( 2 °) 

a=l a=l o=l a=l 

n— 1 n n— 1 n 

/? = r -r • M t =Y^r -r 4- T ia ■ M ia T ia (21) 

a=l a=l a=l a=l 

and where we henceforth suppress the c superscript on vectors and tensors referred to the 
CBF. We have also used I for the identity matrix and M a = R T M*R and M£* = S* T M^S* 
for the moment of inertia of body a in the CBF and IBF a respectively. In the above 
equations, we present two expressions for each of the quantities M, a M , and h^ v . The first 
expression involves quantities referred entirely to the CBF. The second expression is slightly 
more complex but has the advantage that all of the dependence on S a (q) and r a (q) is shown 
explicitly by S a , r a , r a ^, and t™. Note that is the moment of inertia of body a in its 
own IBF and as such is a constant matrix independent of both the shape and orientation. 
By rearranging the terms in Eq. (|18|), we put the kinetic energy in the form 

T= 1 -(u> + A^) • M(« + A v q v ) + "S^T, (22) 

where 



= M"^, (23) 
9^ = K„ - ■ M _1 a I/ . (24) 



Converting the velocities to momenta, we find 



dT 

duj 

dT 



J = — = M(w + A^), (25) 



The vector J is the total angular momentum of the complex. It satisfies the usual body- 
referred "anomalous" commutation relations {Ji, Jj} = — ^^^ijkJki "where is the usual 
Levi-Civita symbol. The classical kinetic energy is now expressible in terms of momenta as 

T = Ij ■ M- X J + I(p„ - A M ■ J)<r(^ - A u ■ J), (27) 
where is the inverse of g^ u . 
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Equation fl27|) decomposes the kinetic energy into two terms. The first term is the kinetic 
energy the complex would have if it were a rigid body of fixed shape. We regard this term 
as the rotational kinetic energy of the complex. The second term we regard as the internal, 
or vibrational, kinetic energy. Often in such rotation-vibration decompositions, a different 
rotational term appears which contains a modified moment of inertia tensor and a modified 
angular momentum vector. We will relate the above decomposition to such alternative 
decompositions in the next section. For now, however, notice that the appearance of A M in 
the internal kinetic energy couples the internal degrees of freedom to the angular momentum. 
For this reason, we call A M the Coriolis potential. Furthermore, we call g^ v the internal metric 
because it acts to square — A M • J. Note that the internal kinetic energy has the same 
"|p — eA| 2 " form as the kinetic energy of a particle in a magnetic field, where the role of 
the vector, or gauge, potential is played by the Coriolis potential and the role of the electric 
charge is played by the angular momentum. 



The kinetic energy given in Ref. |13| for a collection of point particles has exactly the 



same form as Eq. fl27|). However, for collections of point particles, the quantities defined in 
Eqs. flT9| ) - ( f2T| ) do not contain the terms with The appearance of these terms and, 

of course, the introduction of an extra 3n internal coordinates are the sole modifications 
necessary to augment the kinetic energy of a system of point particles to include the rotational 
kinetic energy of the monomers. The fundamental reason why the form of the kinetic energy 
is the same for these two cases is the rotational symmetry of the kinetic energy operator; 
Eq. fl27j ) is in fact a general result, valid for any SO (3) invariant metric. We will discuss this 
matter further in a future publication. 

We comment now on how the preceding results are generalised to include collinear 
monomers and point particles. First, the number of coordinates changes. A noncollinear 
body requires three Euler angles to fully specify its orientation, whereas a collinear body 
requires only two spherical coordinates to specify its orientation (the direction of its collinear 
axis). A point particle, of course, requires no orientational coordinates. Therefore, instead 
of 6n — 3 coordinates as before, there are 3n + 3n n + 2n c — 3 coordinates parametrising the 
centre of mass system. Here, n = n n + n c + n p is the total number of monomers, where n n 
is the number of noncollinear monomers, n c is the number of collinear monomers, and n p is 
the number of point particles. 

We next describe the form of the rovibrational kinetic energy when incorporating collinear 
bodies and point particles. First, the computation leading from Eq. (||) to Eq. (|27D is 
essentially unchanged by the inclusion of collinear bodies and point particles, so long as 
one takes the moment of inertia tensor M a of a point particle to be 0. The moment of 
inertia tensor of a collinear body is explicitly M a = n a (\ — n a n^), where n a is a unit 
vector pointing along the collinear axis and n a is the single nonzero principal moment. The 
quantity n a is a more natural measure of the orientation of a collinear body than S a , since 
S a overparametrises the orientations. We therefore rewrite Eqs. (p~9[ ) - ([H]) in the following 
form, more appropriate for a general complex, 



n-1 



M 



(28) 



a=l 



a=l 



a=n n +l 
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n-1 n n rin+ric 

a M = J^r a x r a ^ + ^S a M^Tj, + K a ii a x n QjAt , (29) 

q=1 o=l a=n„+l 

n—1 n n n n +n c 

a=l a=l a=n n +l 

where we have ordered the monomers with the noncollinear bodies first, the collinear bodies 
second, and the point particles last. We omit the straightforward proof of these equations 
relying instead on the following two observations. First, the orientational contribution from 
the point particles has dropped out since M a is for such particles. Next, the orientational 
contribution from a collinear body is identical to the contribution of a single Jacobi vector. 
This fact is easily understood by modeling a collinear body as two point particles connected 
by a Jacobi vector. The above equations may be combined with Eqs. (|23|) and ( ^4|) to 
obtain the kinetic energy Eq. ( |2T| ) of a general molecular complex. The only note of caution 
occurs if the entire complex should become collinear, in which case M is not invertible and 
singularities may arise. 

2.2 An alternative form of the kinetic energy 

We rearrange the kinetic energy Eq. (|27| ) to place it in a form which is more common in the 
literature on rovibrational Hamiltonians. We define the modified moment of inertia tensor 
M by 

M = M - h^a^, (31) 

where h^ v is the inverse matrix of h^, and the vector K, often called the angular momentum 
of vibration, by 

K = h^&uPu. (32) 
We use these definitions to place the kinetic energy in the form 

T = i(J - K) • M- X (J - K) + l -p^ v Pv . (33) 

The above equation may be verified by comparing the terms of order J 2 , J , and J° in 
Eqs. ( |3"3"| ) and (P7|). The equality of the respective terms is readily apparent from the following 
identities 

g' MU M~ 1 a u = h^M- 1 ^, (34) 
gilv = + h^{a a -M- 1 a T )h TU , (35) 
M- 1 = M^ + ^A^, (36) 

which follow from Eqs. (|3l|) and (p4[) . 

Equation (|3^) provides an alternative rovibrational decomposition of the kinetic energy, 
in a form common in the literature for rovibrational Hamiltonians. For example, the Wilson- 
Howard- Watson molecular Hamiltonian p8|, p9| is expressed in this manner using the Eckart 
conventions, for which various simplifications occur. 



8 



2.3 Changing the internal coordinates and the collective body 
frame 



When the conventions for the internal coordinates or the CBF are changed, the quantities 
defined in this paper do not, in general, remain invariant. Instead, they transform via a 
precise set of rules. Although the analysis of these rules is not critical to the logical flow of 
this paper, we include an account of them for the following two reasons. An understanding 
of transformation rules facilitates conversion between different sets of conventions. This is 
important in actual problems, where it is not uncommon to utilise more than one coordinate 
or frame convention, especially for large amplitude motions. Also, knowledge of transfor- 
mation properties leads to the definition and consideration of quantities which transform 
in a simple manner (that is invariantly or covariantly). Such quantities often have special 



geometric or physical significance. The review of Littlejohn and Reinsch [ 13 1 contains an 
in depth discussion of transformation properties specialised for systems of point particles. 
Since these results are essentially unchanged by the generalisation to include rigid bodies, 



we simply summarise here the key results from Ref. |L3| and comment on how they are to 
be extended. 

We first introduce the important concept of (/-tensors. We consider a new set of coordi- 
nates q" 1 which are functions of the old coordinates g M . A g-tensor transforms under such a 
change in coordinates by contracting each lower Greek index /i with dq^/dq'" and each upper 
Greek index with dq' v jdq^ . The rank of the tensor is the total number of such indices, 
whether upper or lower. For example, T afJi is a rank one g-tensor because it transforms via 

' - d( f (Vf \ 
a ^ dq'^ TaU) 

where r' a ^ is computed from Eq. (|T5|) using the new coordinates q" 1 and T afl is computed 
using the old coordinates g M . Other rank one g-tensors are r^, a M , A M , and g M . Rank 
zero g-tensors, also called g-scalars, are invariant under coordinate transformations. They 
include S*, S a , R, r^,, , uj s al lj s , M s a , M s , M s , T, J s , and K s , as well as these same quantities 
referred (where appropriate) to the CBF or IBF a. Rank two g-tensors include h^, h 11 ", 
g^, and g^. 

Next, we define the concept of an R-tensor, which is important when changing the CBF. 
We consider a new CBF such that the orientation matrix R' with respect to the new frame 
is related to the old orientation R by 

R' = RU(g), (38) 

where U(g) G 50(3) is a smooth function of g. The coordinates g M are held fixed. The 
quantities r a and S a transform via 



r'a = U T r Q , (39) 
S' Q = U T S a , (40) 

where we have omitted the g dependence. We call r a a rank one R-tensor because it has one 
Latin index which transforms with one copy of U T . In general, an R-tensor transforms by 
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contracting each Latin index with U T as in Eq. fl39|) . The rank of the R-tensor is the number 
of Latin indices it possesses. Other rank one R-tensors include n^, u> a , and J. Rank zero 
R-tensors, also called R-scalars, do not depend on the choice of CBF and include T and q^. 
Rank two R-tensors include M and M a . 

The quantity A M is not an R-tensor. Instead, it has a more complicated transformation 
property 



A', = U A M — 7 (41) 



where 7 x = U T U „. An interpretation of this transformation property as the gauge trans- 



formation law of a non-Abelian gauge potential is given in Ref. We simply note that 
such an analysis motivates the introduction of the Coriolis field strength, 

B^jy = A^^ A^y A^j x Ay. (42) 

The Coriolis field strength is a rank two g-tensor and a rank one R-tensor. We will only use 
the Coriolis field strength briefly in Sect. ^. However, it plays a central role in the gauge 
theoretic approach to rovibrational coupling. 

One should be aware that other quantities which have been introduced are also not R- 
tensors. These include h^ u and M. However, it should be noted that both of these quantities 
have counterparts, g^ v and M respectively, which are R-tensors. For further insight into the 



relationship between these two pairs of quantities, see Ref. |13|| . Various other quantities 
which are not R-tensors include r aAt , a M , K, uj, and p M . 

An important observation is that the two decompositions of the kinetic energy, Eqs. ( f27f ) 
and ([33|), differ in their transformation properties. Specifically, Eq. ( p7|) decomposes the 
kinetic energy into two terms which are q- and R-scalars. Equation (p3|), on the other hand, 
decomposes the kinetic energy into two terms which are g-scalars but not R-scalars. Thus, 
the latter decomposition is dependent on the choice of CBF, whereas the former is not. For 
this reason, we view Eq. fl27|) as the fundamental rovibrational decomposition of the kinetic 
energy. However, depending on the CBF convention, the decomposition of Eq. (|33| ) may 
very well be easier to compute. (See Sect. |[) We refer to Ref. |13[] for further discussion of 
these decompositions. 



3 General expressions for the quantum kinetic energy 
3.1 The unsealed kinetic energy 

We temporarily abandon the specific system of a molecular complex in order to present 
expressions for the quantum kinetic energy of a general system. The results of this section 
are largely similar to previous work of Nauts and Chapuisat |30| , which in turn relies on 



several earlier references. We also note that Van der Avoird et al. have employed a similar 
formalism for the water trimer J7j. Here we will summarise relevant aspects of these results 
to fix notation and to lay the foundation for Sect. f|. Also, we apply the formalism to the 



simple example of a single rigid body, which will be of future use. In Sect. |3.2| we present a 
new approach to scaling the wave function. 
We denote the classical kinetic energy by 
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T = ^n a G ab (x)7r b , (43) 

where x stands for a collection of generalised position variables x a , a = 1, . . . , d (d is the 
number of degrees of freedom), n a , a = 1, . . . , d, are generalised momenta, and G ab are the 
components of the inverse metric tensor G" 1 . When repeated, the indices a, b, c, . . . are 
assumed to be summed from 1 to d, both in Eq. fllBp and the subsequent development. The 
momenta n a are linear combinations of the canonical momenta, 

7r a = C a \x)p b , (44) 

where p a is the momentum canonically conjugate to x a and the C b are components of the 
change of basis matrix. Nauts and Chapuisat |3(| call the more general momenta rt a quasi- 
momenta and reserve the term momenta for what we call the canonical momenta p a . The 
kinetic energy Eq. (|43|) is expressed in terms of the canonical momenta p a by 

T = \vcC a c G ab C b d p d = \p a G ab p b) (45) 

where G ab are the components of the inverse metric G _1 = C T G _1 C with respect to the 
canonical momenta. From now on, we omit the explicit x dependence. 



The quantum kinetic energy T is often expressed in the Podolsky form [[H 

1 1 



T=^-^ 7 =p a VGG ab p b , (46) 
2 VG 

where G = det G and where p a is the momentum operator conjugate to x a , 

«. = < 47 > 

setting h — 1. The operator T may also be expressed using the adjoints of the momentum 
operators, 



or more generally, 



blG ab p b , (48) 



f = U\G ab ^ (49) 
where Tx a is the operator corresponding to the classical momentum 7r a , 

K = C a %. (50) 

As observed by van der Avoird et al. 0, even though the adjoint form of the kinetic 
energy represents the same differential operator as the Podolsky form, Eqs. ([i"S|) and ( |4"9"D 
are convenient for evaluating matrix elements since the adjoint of the momentum operator 
effectively acts on the bra to the left. 
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Equations ( fig) and (f49|) are straightforward consequences of the definition of the adjoint. 
For an arbitrary operator A, the matrix element of with respect to wave functions $ and 
$' is 



where the inner product is defined via 



A$ 1$' 



and the volume element is 



(if 



Gdx 1 ...dx° 



G 



detCI 



cfe 1 . . . dx a 



(51) 



(52) 



(53) 



where G = det G. From this definition, one finds that the momenta p a are not in general 
Hermitian but rather satisfy the following relation 



1 



=p a VG. 
G 



More generally, the momenta n a satisfy 



7r; 



K + 



Pb 



GC n 



(54) 



(55) 



where the square bracket notation indicates that pb acts only on the terms inside the brackets. 
Equations (fifiD, (0), and fl5"0P combine to prove Eqs. (^SJ) and (|4*D|). 

An important illustration of the preceding formalism and one which we shall need later 
is that of a single rigid body. We define Euler angles [x 1 , x 2 , x 3 ] = [a, /3, 7] in the usual 
way by R(a,/3, 7) = R z (a)Rj / (/5)R z (7), where R rotates the space frame into the body 
frame and Rj is a rotation about the ith space axis. The body referred angular mo- 
menta [tti, 7r 2 , 7t 3 ] = [Ji, J 2 , J3] are noncanonical momenta related to the canonical momenta 



\Pi,P2,Ps] = [Pa,Pi3,P^ via |2 



" Jl ' 




Pa 


■h 


= c 


PP 


■h 




. Pi . 



(56) 



where 



cos 7 
sin /3 
sin 7 
sin/3 





sin 7 cos 7 cot /3 
cos 7 — sin 7 cot /3 
1 



(57) 



The classical kinetic energy in terms of the angular momenta is T = J • M _1 J/2 where M = G 
is the body referred moment of inertia tensor, which is independent of the Euler angles. The 
volume element is readily computed from Eq. to be 
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dv = Vdet M sin (5 da d(3 d'y = 87rVdet M dR, 



(58) 



where dR = sin/3 da df3 d'y/(8n 2 ) is the normalised Haar measure on SO(3). The quantum 
kinetic energy is expressed in terms of the operators Jj 



" Jl ' 




—id /da 


■h 


= c 


-td/d/3 


. h . 




—id /d^ 



(59) 



Using the volume element dv, one may verify that Jj is Hermitian. (On a deeper level, J, 
is Hermitian because it is a symmetry of the kinetic energy) Thus, the quantum kinetic 
energy Eq. (f49|) acquires the familiar form 



f = -J ■ M -1 J. 



(60) 



3.2 The scaled kinetic energy 

Often it is useful to multiply the original wave function $ by some real positive function 
S(x) to form a new wave function 



(61) 



Such a scaling produces a new kinetic energy operator acting on the new wave function 
In this section, we derive the form of this new kinetic energy operator. Similar discussions 
are given by Nauts and Chapuisat |K| and Chapuisat, Belafhal, and Nauts [33|. The most 
notable distinction between our approach and these earlier accounts is our introduction of a 
new adjoint, shown in Eq. (|64l). This adjoint allows for a different form for the scaled kinetic 
energy operator shown in Eq. (^) and the associated extrapotential term in Eq. (|68|). 

We note that Eq. flBTD induces a new inner product on the scaled wave functions. We 
denote this new inner product with an S subscript and define it via, 



S 



S 



S 2 



S 2 



-dx 1 ...dx d 



(62) 



Thus, V G/S 2 dx l . . .dx d is the volume element associated with the scaled wave functions. 
The operator adjoint taken with respect to this new inner product will in general be different 
from the adjoint taken with respect to the old inner product. To avoid confusion we will 
denote the new adjoint by A^ s \ These two adjoints are related by the following computation 



it(5)^ = = / , Sil/ 



which summarises as 



S 2 



s 2 



(63) 
(64) 
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Equation fl64|) combines with Eq. (|55|) to yield 

vrP = 7rt-2[^ln5]. (65) 

The scaling of the wave function transforms the kinetic energy operator into T$ = 
5T (1/5). Combining Eqs. © and (||), we find 

fs = ~ ^ S) S) G ab • (66) 

By more or less straightforward commutation of operators in the above equation and using 
Eqs. (|55|) and (|65"D, we arrive at the main result of this section 

Ts = \^G ab h + V s , (67) 

where 

V s = -^(G a6 [7r a ln<S][7r 6 ln«S] + [7rt( s )G a6 [7r 6 ln<S]]) 

= i(G a6 [7r a ln < S][7r b ln«S]-[7rtG a6 [7r 6 ln<S]]). (68) 

Comparing Eq. (|67|) to the unsealed expression Eq. (|9|), we note that the two operators differ 
by the additional scalar term in Eq. ( |67|) and the different adjoints which are used. Thus, 
scaling the wave function may be used to place the adjoint of the momenta in an alternative, 
perhaps more attractive, form, but only at the expense of introducing an extrapotential term 
into the kinetic energy. 



4 The quantum kinetic energy of a molecular complex 

We quantise the classical kinetic energy Eq. ( |27j ) using Eq. (|49| ) derived in the previous 
section. This approach requires the operators p^, Jj, and their adjoints. Since the classical 
momentum p^ is canonically conjugate to g M , the quantised operator p^ has the usual form 
of Eq. (g7D 

The quantised angular momenta Jj satisfy the standard "anomalous" commutation relations 
[ Jj, Jj] = —i ^2, k tijkJk an d of course commute with all rotationally invariant operators, for 
example, 

[Ji, M] = [4 M] = [4 A„] = [4 a„] = [Ji, g^] = [J u V] = [J<M = 0- (70) 

To compute the volume element of Eq. ([53|) , we require explicit coordinates covering all di- 
rections of configuration space. This means defining three Euler angles [0 , 6 2 , 6 3 } = [a, (3, 7], 
describing the collective orientation R, which complement the d — 3 internal coordinates q 11 . 
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Here, d = 3n + 3n n + 2n c — 3 is the dimension of the centre of mass system. We adopt 
the Euler angle conventions used in Section |3|. The noncanonical momenta 7Tj = jj and 
^1* = Vii ~ J ' A M are expressed in terms of the canonical momenta Pi = —id/d9 l and p^ by 



Ji 




<V o - 




" Pi ' 










. Pv . 



(71) 



where C i 3 are the components of the matrix in Eq. (|57|) and the sum over j is implicit. The 
full d x d matrix in Eq. (|7T| ) corresponds to the matrix C a b in Eq. fl50p. Its determinant is 
equal to the determinant of the upper left block alone, that is det C = — 1/ sin 8. The metric 
G a b with respect to the momenta 7Tj = Ji and tt^ = p^ — J • A M is seen from Eq. ( ^7] ) to have 
determinant 

G = c/detM, (72) 
where ^ = detp^. The volume element is therefore 



dv = — — — da d(3 dj dq 1 . . . dq d ~ 3 = J a det M sin 8 da dd d'y dq 1 . . . dq d ~ 3 
det C 



= 8vr 2 \/g det M dR dq 1 . . . dq d ' 3 . (73) 
An identity we will use later is the following alternative expression for G 

G = /idetM, (74) 

where h is the determinant of h^ v . This identity follows from the fact that the change of 
basis connecting Eq. fl2"7|) with Eq. ( |3"3"D is orthogonal. 

Since det M and g are rotationally invariant, their presence in dv is irrelevant for the 
computation of jj. Therefore, the computation of j\ reduces to the case of a single rigid 
rotor examined in Section |3| from which we recall that jj = Ji. Therefore, from Eqs. (p7[) 
and fl49D , we find 

f = I j . M^J + l -{p\ - A, ■ J)g^(p v - A, • J). (75) 

The ordering of the operators J{ with respect to the other factors is irrelevant, on account 
of Eq. (|70D . The ordering of the p^ with respect to g^ v and A^, however, is essential. Note 
that the p^ are not in general Hermitian but rather satisfy 

^^^VPOT^^' <T6) 

as easily seen from Eq. ([54]) and the fact that G = G/(det C) 2 = G(sin 3) 2 . 
We now scale the wave function by a factor 



S = V8^G 1/4 = v / 8^(det M) 1/4 ^ 1/4 (77) 

to obtain a new form of the kinetic energy. First, we note that the transformed volume 
element is 
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dv 



dR dq 1 ... dq 



d-3 



(78) 



The angular momenta Jj are still Hermitian with respect to this new volume element, that 
is jj = Ji, as may be noted from Eq. (^) and the fact that S is rotationally invariant. 
However, since the new volume element contains no q^ dependence in the Jacobian prefactor, 
we have the added benefit that is now Hermitian, that is 



TV 



Therefore, the transformed kinetic energy of Eq. 

U 1 



takes the simple form 



= -3 ■ M- X J + -(p„ - A M ■ J)<T (p„ - A u ■ J) + V s , 



(79) 



10) 



where the extrapotential term may be reduced to 



Vs 



i G -l/4 

2 



r) r) 



1) 



dqt 1 dq- 

We observe that Vs is an R-scalar, but not a g-scalar. Therefore, Vs depends on the choice 
of internal coordinates, but not on the choice of CBF. Further discussion of this matter is 
given in Ref. [T3|. 

The quantum kinetic energy may also be placed in a form analogous to Eq. (|33|). The 
unsealed kinetic energy Eq. ([75]) becomes 



f = ~(J - Kt) • M-^J - K) + \v\W v p v , 



(82) 



where 



K = h^arfu. (83) 
Similarly, the scaled kinetic energy becomes 

f 5 = i(J - Kt(<s)) ■ M-\3 — K) + lf>phrp v + V s . (84) 



5 Example: a monomer- atom complex 

We compute the rovibrational kinetic energy explicitly for a system containing a single 
noncollinear rigid monomer with moment of inertia Mx and a single atom, for example, 
Ar-NH 3 . The kinetic energy of such systems has already been studied by Brocks and van 
Koevenjl]], van der Avoird 0, and Makarewicz and Bauder 0. Our presentation is mainly 
designed to illustrate the formalism of the preceding sections, although we believe that the 
derivation of the Coriolis potential A* , Coriolis field strength B l , and internal metric 
is new. 

We define the CBF by fixing it to the rigid monomer. This implies that the matrix S, 
defining the orientation of the monomer's IBF in the CBF, is constant. We take this constant 
to be the identity, 
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S(g) = I. (85) 

Since there is only one rigid body and one Jacobi vector, we drop all "a" subscripts, except 
on Mi, where the "1" serves to distinguish the moment of inertia of the monomer from the 
total moment of inertia M of the complex. The Jacobi vector r locates the atom with respect 
to the monomer and its components may therefore be chosen as the internal coordinates, 
that is 



1 q 2 ,q 3 } = [r!,r 2 ,r 3 ]. (86) 



Thus, 



Ti,^ = <V (87) 
Furthermore, from Eqs. fllBD and ([35] ), we have 



T M = 0. 



Inserting Eqs. (0), and (|88| ) into Eqs. (^) - (pQ|) , we readily obtain 



M = r 2 l - rr T + Mi, (89) 



( rxr J»)i = E^> r i' ( 90 ) 
3 

where the z superscript on Mi indicates that it is referred to the monomer's IBF (which 
agrees here with the CBF) and is hence a constant matrix. 

We proceed to first construct the kinetic energy Eq. (^) which results here in a simpler 
form than Eq. (p7j). Using Eqs. ( pT| ) and (|32| ) we compute 



M = Mi, (92) 

K = r x p, (93) 

where p = [p\ ) P2)Pz\- These are particularly simple results and together with Eq. (0) yield 
the classical kinetic energy 

T = 1(J - r x p) ■ (Mi)- 1 ^ - r x p) + ip • p. (94) 
The quantum kinetic energy requires the further result 

G = /idetM = detMi, (95) 



which follows from Eq. ( 74]) and shows that G is constant. Hence, from Eq. ([FfJ) it is clear 
that p M is Hermitian with respect to the original inner product. Since K = r x p we find 
that K is also Hermitian with respect to both the original and the scaled inner products. 
Furthermore, the extrapotential term of Eq. (RTF) arising in the scaled kinetic energy vanishes. 
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Thus, both the original and the scaled quantum kinetic energies are identical and each is 



formed by simply replacing p and J in Eq. (EH) with p and J respectively. Our results agree 
with earlier derivations by Brocks and van Koeven [|l| and van der Avoird@. 

To simplify the algebra in constructing the kinetic energy Eq. Q27|), we assume the rigid 
body is a spherical top with M\ = k\. The total moment of inertia tensor given in Eq. ( j89|) 
may be explicitly inverted, with the form 



h/T 1 



-1 + 



-rr 



(96) 



r 2 + k n(r 2 + k) 

and combined with Eqs. fl23f) and ( ^4|) to yield explicit forms for the Coriolis potential and 
the internal metric, 



4' 

A* 



The inverse of the internal metric is 



y— 



r z + k 



e ijfi r j j 



r 2 + k 



(97) 
(98) 



■[(r 2 + k)5^ v - r^r u ]. 



(99) 



Eqs. (|96|) , (|9~7|) , and ( p9|) combine with Eq. (|27| ) to yield an explicit form for the classical 
kinetic energy. As earlier, the quantum kinetic energy, both original and scaled, is obtained 
by simply replacing p and J by their operator counterparts, without the need for Hermitian 
conjugates or an extrapotential term. 

It is interesting to compute the Coriolis field strength defined by Eq. (52), 



( r 2 + K )2 



{r-ir a + 2KS itT ). 



(100) 



As the separation r of the atom from the monomer goes to infinity, the Coriolis field strength 
tends toward B l — > ^^uafira/r^. We change the CBF, as in Eq. ( j38|) , via a matrix U(r) 
which rotates z into r/r. Then, since is a rank one R-tensor, as r goes to infinity, the 
new field strength tensor approaches 



b: 



U T B 



(101) 



The above asymptotic form is that of a (non-Abelian) monopole field [34]. A similar 
monopole field is already known to exist in the three-body problem |T5 , IB , a fact which 
has led to several useful applications [|1^, 0, |1^, . We remark that the above asymptotic 



form is valid even if the monomer is an asymmetric top. 



6 Conclusions 



We have computed the kinetic energy of an arbitrary molecular complex for arbitrary co- 
ordinate and body frame conventions. In so doing, we have tried to provide an efficient 
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framework in which explicit Hamiltonians may be readily computed for specific choices of 
coordinates and frames. We have provided a discussion of transformation properties to facil- 
itate the changing of these conventions. Our formalism is illustrated with the example of a 
monomer-atom system, and more complex systems may be handled with similar ease within 
our framework. 

One of the more novel and intriguing aspects of our derivation is the appearance of the 
Coriolis potential and the various insights which are possible by adopting a gauge theoretical 
viewpoint. We briefly cite two areas of current research which are based on this perspective. 
First, using gauge theoretic reasoning we have managed to generalise the Eckart conditions, 
so often employed for small vibrations in molecules, to systems of rigid bodies. Much of 
the formalism for small vibrational analysis in molecules can then be readily ported over 
to study small amplitude vibrations in clusters of rigid molecules. Second, we have been 
able to understand rotational splittings in molecules with internal rotors as a sort of Coriolis 
Aharonov-Bohm effect. These applications and others will be the subject of future publica- 
tions. 
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